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Lewis  B.  Wetzel 


ON  THE  THEORY  OF  ELECTROMAGNETIC  SCATTERING 
FROM  A  RAINDROP  SPLASH 


I.  INTRODUCTION 

Interest  in  electromagnetic  scattering  by  splashes  derives  primarily  from  interest  in  the  effect  of 
rain  on  radar  sea  backscatter.  Radar  operators  will  tell  you  that  rain  tends  to  suppress  sea  back- 
scatter,  and  the  observations  of  operators  should  always  be  taken  seriously.  Yet  there  has  been  little 
in  the  way  of  reliable,  quantitative  experimental  information  about  the  interaction  between  rain  and 
rain-driven  sea  backscatter,  and  no  real  theory  whatsoever.  Laboratory  measurements  by  Moore,  et 
al,  [1979|  with  artificial  “rain”  suggested  that  for  light  “winds”  the  backscatter  level  increased  with 
rain  rate,  while  for  heavy  winds  rain  made  little  difference.  In  measurements  in  natural  rain  over  the 
Chesapeake  Bay,  Hansen  1 19841  found  that  even  a  light  rain  (2  mm/hr)  changes  the  spectral  character 
of  sea  backscatter  at  moderate  wind  speeds  (6  m  'sec)  by  introducing  a  significant  high  frequency  com¬ 
ponent.  He  also  found  some  evidence  in  support  of  the  radar  operators  at  low  grazing  angles  with 
horizontal  polarization,  which  is  the  most  common  operating  regime  for  Navy  radars.  Since  these 
two  studies  constitute  the  total  published  effort  on  this  problem,  an  understanding  of  the  complex 
phenomenology  of  rain/sea-scatter  interactions  will  have  to  be  approached  essentially  from  scratch. 

This  paper  takes  a  first  step  by  developing  a  model  for  scattering  from  an  isolated  splash  and 
applying  this  model  to  a  statistical  description  of  the  radar  cross  section  of  natural  rain  splashes  on  a 
“calm"  water  surface. 

II.  THE  SPLASH  PROCESS 

Not  everyone  likes  the  rain,  but  there  are  few  who  have  not  wondered  at  the  unexpected  display 
of  a  raindrop  falling  on  a  water  surface  —  the  explosive  jet  rising  from  the  center,  the  fleeting  silver 
ring  that  vanishes  almost  before  you  realize  what  it  is.  Over  100  years  ago,  his  fascination  with  this 
event  led  Worthington  1 1 882.  1962|  to  develop  a  remarkable  photographic  system  for  recording  the 
process  with  millisecond  precision,  using  “state-of-the-art”  equipment  of  his  day.  His  raindrop  was  a 
ball  of  water  rolled  off  a  smoked  watchglass.  activated  by  a  magnetic  latch  which  also  opened  the 
shutter  of  a  wet-plate  camera  whose  emulsion  he  had  freshly  mixed.  After  a  suitable  delay,  a  bank  of 
Leyden  Jars  was  discharged  into  an  arc  to  illuminate  the  scene  with  a  flash  of  millisecond  duration. 
Some  examples  of  his  wizardry  arc  seen  in  Fig.  I,  where  we  see  the  three  basic  structures  associated 
with  a  splash:  the  “crown,”  the  “stalk"  (or  “plume"  or  “jet”),  and  the  “ring  wave.”  Since  then, 
of  course,  the  beauty  and  novelty  of  the  event,  captured  easily  by  modern  high-speed  flash  photogra¬ 
phy.  have  made  stop-action  and  slow— motion  photos  of  splashes  familiar  to  almost  everyone.  A  gal¬ 
lery  of  photos  taken  recently  by  Cavaleri  [1985]  is  shown  in  Figs.  2  and  2. 

It  is  clear  from  the  photographic  evidence  that  the  splash  process  produces  three  well-defined 
structures  whose  dimensions  seem  to  be  related  to  the  drop  diameter  I)  They  occur  in  time 
sequence,  with  the  entire  process  proceding  as  follows:  First  a  “crown”  with  an  irregular  upper 
boundary  grows  from  the  initial  impact  to  a  height  of  about  I  I)  and  diameter  of  3-4  1).  The  growth 
and  decay  of  the  “crown"  lasts  about  50  milliseconds,  leaving  an  elevated  ring  of  water  moving  out 
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from  the  site  with  a  central  depression  at  its  center.  The  floor  of  this  central  depression  rises  into  a 
kind  of  cone,  from  which  grows  a  “stalk”  carrying  a  ball  of  liquid  at  its  top.  The  stalk  diameter 
seems  to  be  about  3/4  D  and  reaches  a  maximum  height  of  4-6  D.  Its  growth  and  decay  lasts  100- 
200  msec,  and  as  it  reaches  its  peak  and  starts  to  decay,  the  ball  may  be  squeezed  off,  with  the 
column  separating  into  several  droplets  (Fig.  3)  —  or  it  may  not  (Fig.  2).  The  collapsing  stalk  pro¬ 
duces  a  second  annular  wave  that  propagates  outward  in  pursuit  of  the  original  “crown”  wave,  the 
whole  thing  preceded  by  what  appears  to  be  a  set  of  waves  of  very  small  wavelength  (parasitic  capil¬ 
laries?)  If  secondary  droplets  have  been  squeezed  off  in  the  stalk,  they  fall  back  to  form  much 
weaker  systems  of  secondary  “ring  waves.”  Variations  in  the  details  of  the  splash  process  are  prob¬ 
ably  due  to  variations  in  such  parameters  as  surface  and  drop  temperatures,  surface  cleanliness,  drop 
size  and  velocity.  For  large,  high-speed  drops,  for  example,  the  walls  of  the  “crown”  are  seen  to 
meet  over  the  top  of  the  initial  crater  to  form  a  “bubble.” 

The  detailed  hydrodynamics  of  the  splash  process  is  not  well  understood.  In  fact,  the  major 
reference  for  rain-like  splashes  remains  the  mostly  descriptive  work  by  Worthington  cited  above.  The 
more  recent  work  on  splashes  is  directed  toward  special  cases  —  e.g..  high-velocity  impacts  (Kngel, 
1966|.  and  low-angle  sprays  in  which  drops  bounce  off  the  surface  like  skipping  stones  |Jayaratne  and 
Mason,  1964J.  However,  some  insights  into  general  splash  dynamics  may  be  found  in  these  refer¬ 
ences.  For  example,  in  Engel’s  experiments  the  falling  drop  was  dyed  red,  and  white  particles  were 
suspended  in  the  target  water  to  help  visualize  the  resulting  internal  flow  patterns.  It  was  found  that 
the  ball  atop  the  “stalk”  contained  the  red  dye  of  the  original  falling  drop,  and  that  the  growth  and 
decay  of  the  crown  and  the  emergence  of  the  stalk  were  associated  with  an  oscillating  toroidal  vor- 
ticity.  But  knowledge  of  the  physics  of  splashes  is  really  not  necessary  for  the  development  of  scatter¬ 
ing  models.  It  is  sufficient  to  to  know  the  structural  forms  of  the  various  parts  of  the  splash,  and 
photographic  evidence  of  the  type  seen  in  the  Figures  above  can  provide  this  information. 

III.  SPLASH  SCATTERING  EXPERIMENTS 

In  the  first  recorded  measurement  of  its  kind,  Hansen  |1984|  measured  radar  backseatter  versus 
rain  rate  from  a  “calm”  natural  water  surface.  This  was  done  over  the  Chesapeake  Bay  using  an  X- 
band  radar  with  both  horizontal  and  vertical  polarizations.  A  patch  of  calm  water  surface  was  under 
illumination  just  as  it  started  to  rain,  and  as  the  rain  intensity  steadily  increased,  the  record  shown  in 
Fig.  4  was  obtained.  While  it  is  dangerous  to  generalize  from  a  single  measurement,  it  is  clear  that 
for  this  rain  sample,  at  least,  there  was  a  strong  polarization  dependence  at  low  rain  rates,  while  at 
the  higher  rain  rates,  both  polarizations  reached  backseatter  levels  equal  to  those  for  a  moderate 
wind-driven  sea  at  the  same  grazing  angle  (about  -40  dBnr/nr  at  2  grazing.) 

In  order  to  determine  the  scattering  cross  sections  of  individual  splashes,  Hansen  set  up  a 
laboratory  experiment  in  which  scattering  from  splashes  produced  in  a  laboratory  tank  by  drops  of 
known  size  was  measured  using  a  high-resolution  time  —  domain  re  Hectometer  |I985.  1986)  This 
instrument  recorded  the  time  history  of  scattering  from  the  entire  splash  process  on  a  split-screen 
display  ,  which  showed  this  history  updated  at  17  msec  intervals  along  with  a  snapshot  of  the  instan¬ 
taneous  state  ot  the  developing  splash.  An  example  of  such  a  sequence  is  given  in  Fig.  5  for  verb 
cally  polarized  backseatter.  Figure  6  shows  two  pairs  of  records  of  the  complete  scattering  process, 
with  V  Pol  on  the  top  and  H  Pol  below,  the  first  for  a  drop  of  4mm  diameter  and  the  second  for  a 
3mm  drop.  Drops  of  the  same  size  tended  to  give  similar  results  (observe  the  close  similarity 
between  the  scattering  profiles  in  Figs.  5  and  6a  for  the  same  size  drop),  so  the  differences  between 
the  4mm  and  3mm  returns  indicate  the  nature  of  drop-size  sensitivities  in  the  scattering  process  It 
should  be  noted  that  the  two  polarizations  could  not  be  recorded  simultaneously  for  the  same  drop,  so 
the  upper  and  lower  records  in  Fig  6  belong  to  (wo  different  drops  of  the  same  size.  We  will  be  dis 
cussing  these  records  in  greater  detail  in  the  next  section 
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This  meager  bcxiy  of  scattering  data  is  about  all  there  is  at  the  present  time,  and  the  laboratory 
data  cannot  really  be  considered  representative  of  natural  rain,  since  the  drops  fell  through  distances 
short  of  those  reauired  to  reach  terminal  velocity  in  the  open  air.  Nevertheless,  the  splashes  produced 
in  the  laboratory  and  in  nature  look  very  similar,  and  it  is  the  look,  rather  than  the  detailed  physics, 
of  the  splash  process  that  we  will  use  to  guide  our  nuxieling  efforts. 


IV.  MODELING  CROWN  AND  STALK  SCATTER 


Both  crown  and  stalk  have  cylindrical  symmetry,  so  it  is  tempting  to  model  them  as  finite  water 
cylinders  standing  out  of  an  infinite  water  surface.  A  full-scale  boundary  value  problem  even  for  this 
simple  geometry  would  be  quite  difficult  to  solve,  so  we  revert  to  an  even  simpler  ad  hoc  model  in 
which  the  cylinders  are  quasi-metallic  (impenetrable)  and  the  water  surface  enters  simply  as  a  reflect¬ 
ing  plane  with  a  prescribed  relection  coefficient 


The  basic  scattering  geometry  is  shown  in  Fig.  7.  where  the  scattering  cylinder  is  represented  by 
a  stack  of  discs,  each  being  a  slice  of  an  infinite  cylinder  illuminated  front,  and  scattering  into,  the 
pair  of  possible  paths  created  by  the  presence  of  the  surface.  The  scattered  field  for  each  disc  is  writ¬ 
ten  in  terms  of  radial  eigenfunction  expansions  appropriate  to  illumination  of  an  infinite  cylinder  by 
plane  waves  arriving  along  the  direct  and  surface-reflected  paths  and  returning  to  the  source  along 
these  same  paths.  Then,  as  described  in  the  Appendix,  the  total  scattered  field  is  assembled  by 
integrating  over  the  stack  of  discs.  This  approach  is  an  elaboration  of  the  procedure  used  by  Kerr 
[1952|  to  calculate  the  scattering  from  metallic  cylinders  ignoring  end  effects.  Kerr  treats  an  isolated 
perfectly  conducting  finite  cylinder  of  length  h  illuminated  broadside  (i/-  =  0)  by  either  vertically 
(£[  =  £,.<".)  or  horizontally  (E„  =  £//<’,)  polarized  waves.  The  expressions  given  by  Kerr  for  the 
cross  sections  are  well  known: 
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where  J„  is  the  ordinary  Bessel  function,  Hn  is  the  (outgoing  wave)  Hankel  function.  A  =  27r/X.  anil 
6„„  is  the  Kronecker  delta.  By  the  procedure  described  in  the  Appendix,  a  multipath  factor  is  intro 
duced  to  account  tor  the  reflecting  surface,  and  an  additional  factor.  D‘.  accounts  for  the  fact  that  the 
cylinder  is  made  of  water.  With  these  modifications.  (I)  and  (?)  take  the  form: 
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depends  on  the  complex  reflection  coefficients  tor  the  water  surface  and  is  given  in  ( A 5 ) . 


and  ,,  are  empirical  factors.  A  further  assumption  required  in  the  conversion  of  (I)  and  (2)  into 
(3)  and  (4)  is  that  at  the  disc  surfaces,  the  exterior  boundary  conditions  tor  a  highly  reflective  surface 
"resemble  ”  those  for  a  perfectly  conducting  surface.  Solutions  for  a  dielectric  cylinder  [see  King  and 
Wu.  1959.  C'hapt  2[  would  indicate  that  this  is  a  good  assumption  for  thin  water  cylinders  More 
over,  in  comparing  the  backscatter  cross  sections  of  metal  and  water  spheres.  Aden  [King  and  Wu. 
loc.  tit  |  found  that  the  scattering  behavior  versus  ka  was  virtuallv  the  same,  both  theoretically  and 
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experimentally,  only  the  mean  cross  section  level  for  the  water  spheres  was  about  half  that  for  the 
metal  spheres.  This  result  supports  our  '’resemblance”  assumption,  and  suggests  a  value  of  about  0.5 
for  the  empirical  factor  D:.  (The  reflection  coefficient  of  a  plane  water  surface  is  about  0.75  at 
microwave  frequencies,  further  supporting  this  choice.)  The  wall  of  the  "crown”  is  seen  to  be  quite 
thin  in  some  of  the  photographs,  so  we  might  be  concerned  whether  a  solid-post  model  is  suitable  for 
this  structure.  However,  calculations  of  the  reflection  coefficient  for  a  thin  sheet  of  water  |see  Strat¬ 
ton,  1941,  p.  5151  indicate  that  at  X-band,  at  least,  the  reflection  coefficient  remains  quite  high  for 
thicknesses  down  to  about  0.25mm,  which  is  certainly  very  thin. 

The  cross  sections  can  be  expressed  more  compactly  by  letting  (Ah r)  times  the  respective  sum 
terms  in  (3)  and  (4)  be  denoted  by  the  symbol  CV  H,  putting  D:  =  (1/2),  and  writing 

aVM(a,h-f\i)  =  (1/2)  Cv  H  (ka  cos  i)  ■  hl  ■  \  Fv  "  (f .  kh  )  \ '  (5) 

where  the  dependence  on  cylinder  radius,  a,  height,  h,  and  illuminating  frequency,  f.  and  grazing 
angle  i/-,  is  made  explicit.  The  cylindrical  scattering  factors  C1  11  were  computed  from  their  eigen¬ 
function  expansions,  with  the  results  shown  in  Fig.  8.  The  multipath  factor  Fl  11  depends  on  fre¬ 
quency  through  both  kh  and  the  complex  reflection  coefficient  Rv(f  for  horizontal  polarization. 
RH  ~  - 1  for  almost  all  frequencies  at  relatively  low  grazing  angles.  The  factors  |  Fy  11  | :  are  plot 
ted  in  Fig.  9  for  sea  water,  using  the  amplitude  and  phase  of  Rv  at  3  cm.  (X-band)  given  in  Saxton 
and  Lane  1 1955] . 

A.  Application  of  the  Metallic  Cylinder  Model: 

We  will  begin  by  applying  the  simple  scattering  model  outlined  above  to  the  crown  and  stalk  phases 
of  a  splash  event  characterized  by  the  parameters  of  the  particular  measurement  recorded  in  Fig.  6a: 
X-band  radar  (10  Ghz)  viewing  the  splash  of  a  4  mm  drop  at  a  grazing  angle  of  15°.  It  will  be 
necessary  to  make  certain  assumptions  about  the  dimensions  and  time  behavior  of  the  splash  struc¬ 
tures.  Both  crown  and  stalk  will  be  assumed  to  rise  and  fall  linearly  over  their  lifetimes,  the  crown 
lasting  50  msec  and  rising  to  a  maximum  height  of  ID  (4mm)  with  a  mean  diameter  of  3D  (12mm), 
while  the  stalk  lasts  150msec,  rising  to  6D  (24mm)  with  a  diameter  of  (3/4)D  (3mm).  Thus  in  Eqn. 
(5)  the  crown  will  have  ka,.  =  1.3  with  a  maximum  kh,.  =  0.8.  and  the  stalk  ku ,  -  0.3  with  a  max¬ 
imum  kh,  =  5.  The  resulting  time  histories  of  the  crown  and  stalk  cross  sections  are  plotted  in  Fig. 
10  for  both  horizontal  and  vertical  polarizations.  For  scaling,  we  have  shown  the  return  from  a  4mm 
spherical  water  drop  as  it  would  be  seen  at  the  peaks  of  the  interference  pattern  above  the  surface  for 
both  polarizations  (V-Pol  line  at  -14dBcnr;  H-Pol  at  -7dBcnr).  Their  difference  is  due  to  the  dif¬ 
ferent  surface  reflection  coefficients  for  the  two  polarizations. 

In  comparing  the  model  predictions  in  Fig.  10  with  the  actual  measurements  in  Fig.  6.  the  almost 
total  disagreement  for  the  vertically  polarized  stalk  returns  is  most  striking.  But  there  are  also  some 
strong  similarities  that  become  clearer  if  Fig.  10  is  plotted  against  the  same  linear  power  scale  that 
was  used  in  Fig.  6.  This  is  done  in  Fig.  11a,  using  the  peak  of  the  vertical  crown  return  as  a  refer¬ 
ence  for  comparison  with  the  first  200msec  of  the  experimental  curve  sketched  below  it  (Fig.  11b. 
sketched  from  Fig.  6a.)  The  small  returns  on  the  left  in  the  experimental  records  are  ascribed  by 
Hansen  to  the  passage  of  the  bare  drop  through  the  interference  pattern  above  the  surface  before  it 
strikes,  so  these  levels  should  correspond  to  the  "bare  drop”  cross  sections  given  in  Fig.  10.  Bare 
drop  returns  are  simulated  on  Fig.  11a  by  short  vertical  lines  at  the  approximate  locations  one  would 
expect  them  to  occur  relative  to  the  splash  process.  There  is  a  bit  of  ambiguity  in  properly  identify 
mg  the  small  pips  in  this  part  of  the  records  (for  example.  Hansen  has  suggested  lhal  the  right  hand 
pip  on  the  H-pol  record  might  be  the  crown  return),  but  it  they  could,  indeed,  be  validated  as  bare 
drop  returns  the  measurement  would  become  self  calibrating,  since  the  cross  section  of  a  water  sphere 
of  given  diameter  is  rather  precisely  known.  However,  if  we  accept  all  of  the  assumptions  made  thus 


far,  we  would  conclude  from  Fig.  11  that  the  simple  “metallic  cylinder"  model  is  in  surprising 
agreement  (within  a  few  dB)  with  the  measured  behavior  of  Crown  scattering  and  Horizontal  Stalk 
scattering,  but  fails  to  account  for  the  sharp  cut-off  and  oscillating  behavior  of  the  Vertical  Stalk 
return.  Obviously,  for  these  large  splashes,  there  is  something  seriously  wrong  with  the  assumptions 
underlying  our  cylindrical  scattering  model  for  vertical  polarization. 

B.  A  Dielectric  Waveguide  Model  of  the  Stalk: 

In  applying  the  metallic  cylinder  model,  the  basic  assumptions  are  that  the  scattered  field  components 
just  outside  the  surface  of  the  cylinder  are  uniform  in  the  axial  direction,  that  there  are  neither  interior 
fields  nor  end  effects,  and  that  the  aqueous  nature  of  the  cylinder  is  accounted  for  by  the  empirical 
factor  D‘=0.5.  However,  dielectric  cylinders  can  support  a  variety  of  internal  waveguide  modes, 
and  if  these  were  excited  by  an  incident  wave  one  might  expect  some  interesting  scattering  behavior. 
Figure  12a  shows  the  field  configurations  in  a  dielectric  rod  for  the  first  (lowest)  two  waveguide 
modes:  TMm  and  TE0).  The  cut-off  freouenvies  for  both  modes  are  given  by 

/,  =  1  1 .5/M A'  -  Do  GHz.  (6) 

where  K  is  the  dielectric  constant  of  the  rod  and  a  is  its  radius  in  cm.  |see  Jackson,  1962,  or  John¬ 
son,  1965].  For  water  at  X-band  frequencies,  K  =  60,  so  the  cut-off  frequency  for  a  3mm  stalk 

diameter  is  9.9  Ghz.  Although  we  will  not  hold  these  formulas  to  exact  applicability,  it  is  comforting 

that  the  cut-off  frequency  lies  in  the  right  range. 

While  the  TE(M  fields  would  be  very  difficult  to  excite  with  an  external  plane  wave,  the  magnetic  field 

in  the  TM01  mode  is  a  simple  elaboration  of  the  field  in  the  metallic  cylinder  model.  This  is  illus¬ 
trated  in  Fig.  12b  where  the  magnetic  field  circling  the  rod  at  the  right  is  uniform  along  the  rod, 

while  that  belonging  to  the  TM0,  mode  on  the  left  has  the  same  symmetry',  but  changes  direction 

along  the  rod  with  the  periodicity  of  the  wavelength  in  the  rod.  This  wavelength  is  obtained  from  the 
dispersion  diagram  for  the  mode,  shown  in  Fig.  13  (based  on  Johnson.  1965,  Fig.  4.45). 

In  order  to  explore  the  implications  of  this  idea,  we  must  imagine  that  somehow  the  vertically  polar¬ 
ized  incident  wave  excites  a  TM,,,  mode  as  a  standing  wave  in  the  growing  stalk.  The  azimuthal 
magnetic  field  in  this  mode  is  continuous  across  the  surface,  and  enters  the  scattering  integral  in  (Al) 
to  become  the  source  of  the  scattered  field.  We  see  in  Fig.  14a  that  as  the  stalk  grows  from  zero,  the 
polarity  of  the  magnetic  field  changes  with  each  half-wavelength  of  stalk  length.  This  phase  depen¬ 
dence,  when  added  to  that  produced  bv  oblique  incidence  and  the  grazing  angle  dependence  of  the 
surface  reflection  coefficient,  could  give  rise  to  a  scattered  field  of  considerable  complexity.  Without 
attempting  an  actual  solution  of  the  dielectric  rod  problem,  we  will  simulate  the  axial  behavior  of  the 
surface  magnetic  field  by  inserting  a  factor  with  the  right  periodicity  into  the  integrand  of  (Al).  A 
cosine  function  is  chosen  to  ensure  that  the  electric  field  at  the  base  of  the  rod  will  be  tangential  to 
the  rod  surface,  resulting  in  the  substitution: 

H(r')  -  H(r’\  cos  (kvz').  (7) 

where  kv  is  the  wavenumber  in  the  dielectric  rod.  as  given  In  the  dispersion  curve  in  Fig.  13.  This 
cosine  factor  now  appears  under  the  integral  in  (A3),  with  the  integrated  expression  resembling  (A5) 
but  having  twice  as  many  terms  In  order  to  determine  the  value  to  use  for  ky ,  we  return  to  the  cut¬ 
off  relation  (6).  As  noted,  the  cut  off  frequency  for  a  3mm  diameter  water  rod  is  at  X-band.  so  w'e 
expect  the  guide  wavelength  to  differ  very  little  from  the  free  space  wavelength  at  the  measurement 
frequency  Using  the  cosine  from  (7)  in  (A3)  with  k i;  k„.  and  performing  the  integral  to  find  the 
equivalent  ol  (A5)  for  the  dielectric  roil  problem,  we  obtain,  finally,  the  scattering  behavior  plotted  in 
Fig.  14b  for  the  15  grazing  angle  used  in  the  experiment.  The  experimental  profile  is  sketched  to  the 
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same  scale  to  compare  the  time  histories  (quantitatively  it  turns  out  to  be  about  twice  the  predicted 
value.)  Clearly,  the  dielectric  rod  model  looks  promising  as  an  explanation  for  the  peculiar  behavior 
of  the  vertically  polarized  splash. 

Additional  support  for  this  hypothesis  is  found  in  some  measurements  being  made  by  Hansen  (private 
communication)  in  which  the  equipment  used  to  obtain  the  records  shown  in  Figs.  5  and  6  measured 
the  backscattcr  from  a  vertical  water  column  contained  in  a  plastic  tube.  By  recording  the  time  his- 
tory  of  the  scattered  signal  as  the  water  was  drained  rapidly  from  the  tube,  the  scattering  behavior  of 
a  growing  (or  decaying)  •'stalk"  was  simulated.  An  example  is  shown  in  Fig.  14'a  for  a  9cm  column 
with  a  diameter  of  6mm,  illuminated  at  a  grazing  angle  of  27°.  In  Fig.  14'b  the  corresponding 
theoietieal  result,  based  on  using  (7)  in  (A3),  is  plotted  for  kv  -  k:l  While  the  similarity  of  the  two 
curves  is  striking,  there  are  problems  in  both  the  experiment  (a  water  film  is  left  inside  the  plastic 
tub;  during  the  brief  3(X)nisec  draining  time),  and  the  theory  (the  mode  structure  in  the  water  column 
would  be  complicated  by  the  plastic  sleeve  supporting  it.  m  fact,  it  was  necessary  to  use  the  free 
space  wa.eiength  in  the  rod  to  obtain  the  agreement  shown,  even  though  the  guide  wavelength  for  the 
6mm  -•uimm  would  have  been  much  shorter).  Nevertheless,  the  agreement  shown  in  Figs.  14  and 
14'  certainly  enhances  the  plausibility  ol  the  model. 

The  asvmmetries  in  the  scattering  profiles  for  the  falling  drops  are  probably  due.  at  least  in  part,  to  a 
difference  between  the  growth  and  decay  times  for  the  stalk:  the  water  column  must  labor  to  push 
itself  upward  against  the  forces  of  gravity  anil  surface  tension,  and  once  having  used  up  the  energy 
axaiiui'l.  oom  the  falling  raindrop,  it  "crashes"  in  a  much  shorter  interval  (sometimes  breaking  up  in 
fnc  p.oic^.  as  in  Fig.  3).  The  displacement  of  the  central  peak  in  Fig.  14b  and  the  rounded  leading 
edge  m  me  H  pol  return  in  Fig.  6  give  evidence  of  this  type  of  asymmetry.  There  are  several  other 
featuiC'  m  the  experimental  records  that  require  explanation.  The  large  return  at  about  260  msec  in 
the  V  p  i  record  of  Fig  6a  (see  also  Fig.  5)  is  associated  with  a  secondary  stalk  formed  by  the 
"crash''  of  the  main  stalk.  Although  these  secondary  stalks  can  be  quite  tall,  they  are  observed  to  be 
very  mm  t  hus  their  cut-off  frequency  would  be  much  too  high  for  dielectric-rod  scattering,  and  they 
would  probably  scatter  in  the  "metallic  cylinder"  mode.  In  fact,  this  V-pol  return  closely  resembles 
the  metallic  cylinder  H-pol  return  in  the  figure  below  it.  Moreover.  Hansen's  preliminary  measure¬ 
ments  of  scattering  by  water  columns  showed  a  thinner  column  of  the  same  height  as  that  in  Fig. 
14  a  io  scatter  more  strongly  than  the  fatter  column,  and  to  display  greatly  reduced  oscillations  dur¬ 
ing  its  decay . 

The  i!  m  i  equally  spaced  returns  marching  off  to  the  right  are  produced  by  interference  between  the 
adv.tikiiic  and  receding  edges  of  the  "ring  wave."  as  we  will  find  in  the  next  section.  The  3mm 
stalk  ml.  6b  could  be  viewed  as  a  variant  ot  the  4mm  ease,  or  something  quite  different.  We  have 
-.ecu  ••■...  .Hilarities  between  two  different  4mm  drops  in  Figs.  5  and  6a.  but  we  have  no  additional 
data  :  1  ;  'mm  diop 


V-iulc  . -.w  cannot  avoid  leelmg  at  this  point  that  scattering  from  the  crown  and  stalk  phases  of  the 
lam.!:  p  •;  ash  is  tairiv  well  understood,  at  least  qualitatively,  confidence  in  numerical  cross  section 
picon'"  aiII  have  to  wait  until  measurements  ot  splash  scattering  are  available  for  a  wider  range  ot 
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led  earliei  the  "King  Wave"  is  a  fugitive  event,  appearing  a  fraction  of  a  second  alter 
1  earning  aw, iv  a  portion  ot  the  encigv  ot  the  lalleti  raindrop  In  this  sense,  the  splash 
,  he  viewed  as  in  inter  medial  v  lor  convening  the  k  met  is  encrgv  ot  rain  into  the  energy  of 
surer  v\. ives  covciing  the  surlace.  although  a  lough  calculation  shows  that  only  a  small  part 
awav  in  the  line  waves.  I  he  chainclorisiics  ol  ihc  individual  Rime  Waves  mav  best  be 
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determined  by  high-speed  photography,  or  what  is  almost  as  good,  by  using  a  video  camera.  Figure 
15  is  a  sample  of  a  frame-by-frame  video  record  of  natural  rain  falling  on  a  water  surface  |Cavaleri, 
1985-86].  The  ring  waves  arc  seen  to  stabilize  quickly  into  a  simple,  well-defined  shape  that  is 
retained  even  as  the  waves  pass  through  each  other.  This  is  one  of  the  properties  often  ascribed  to  a 
Soliton,  although  any  linear  wave  group  should  also  display  this  interpenetrability. 

A.  Characteristics  of  the  Ring  Wave: 

Although  of  very  low  resolution  compared  to  photographs,  the  video  frames  in  Fig.  15  contain  a 
surprisingly  large  amount  of  information.  The  wave  shape,  and  hence  its  peaks  and  troughs,  can  be 
inferred  from  the  shadowy  relief.  The  maximum  stalk  height  provides  a  rough  yardstick  for  measur¬ 
ing  lengths,  w  hile  the  video  frame  rate  of  one  frame  per  17  msec  is  the  tick  of  the  clock.  Using  these 
measures,  we  should  be  able  to  obtain  estimates  of  the  waveform,  amplitude  and  velocity  of  the  ring 
waves.  Six  sequences  of  the  type  shown  in  Fig.  15  were  analyzed  in  the  manner  of  Fig.  16.  The  ini¬ 
tial  wave  peak  emerges  with  the  collapse  of  the  crown  in  frame  3.  The  stalk  reaches  its  maximum 

height  in  frame  5.  and  as  it  collapses  it  produces  the  second  ring  wave  peak  in  frame  7.  In  this 

sequence  we  see  a  second  stalk  rising  and  collapsing  in  frames  9-13.  giving  rise  to  a  weak  secondary 

ring  wave  starting  in  frame  13.  (We  have  ignored  the  considerably  weaker  parasitic  capillaries  that 

preceed  the  initial  ring  wave  component.)  All  of  the  sequences  analyzed  looked  the  same,  although 
not  all  displayed  the  "second  splash."  The  dashed  lines  trace  the  propagation  of  the  wave  peaks,  and 
their  slopes  give  the  wave  velocity.  Rough  observation  indicated  that  the  maximum  stalk  heights 
were  of  the  order  ot  1  inch,  or  about  2.5cm.  Using  this  value  to  scale  the  radius  of  the  expanding 
ring  wave,  the  velocities  were  of  the  order  of  25  cm/sec,  which  is  close  to  the  minimum  of  the 
velocity/wavelength  characteristic  for  water  waves  (the  transition  between  gravity  and  capillary 
waves.)  The  basic  waveform  produced  by  the  sequential  collapse  of  crown  and  stalk  consists  of  two 
peaks  separated  by  a  trough.  Applying  the  scale  used  to  obtain  the  wave  velocity,  the  two  peaks 
appear  to  be  separated  by  about  1. 5-2.0  cm.  which  lies,  again,  at  the  minumum  of  the 
velocity /wavelength  characteristic  for  water  waves.  This  is  reassuring,  because  it  is  just  the  kind  of 
natural  response  that  one  would  expect  for  an  impulsive  excitation  of  this  type. 

The  waveform  will  be  approximated  by  a  circular  wave  traveling  with  velocity  V  and  having  the  pro¬ 
file  of  a  double-  Gaussian  of  the  form 
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where  p,,  is  a  reference  radius,  and  the  parameters  A.  B,  and  C  are  chosen  to  provide  a  zero-mean 
water  waveform,  which  requires  that  C~A'B  The  cross-sectional  contour  of  such  a  waveform  with 
peaks  separated  by  2  cm  is  shown  in  Fig.  17.  and  seems  to  reflect  rather  well  the  general  impression 
of  the  relatively  deep  trough  separating  the  two  peaks.  For  this  particular  waveform.  A  =  lcm. 
B  65cm.  and  ('  --  1.54:  we  will  use  these  values  in  die  numerical  example  below. 

B.  A  Perturbation  Theory  of  Ring  Wave  Scattering: 

It  is  obvious  from  the  Figures  lhal  the  amplitude  of  the  ring  wave  can  he  at  most  a  millimeter  or  so. 
making  this  one  of  the  few  scattering  features  on  a  water  surface  to  which  a  perturbation  approxima¬ 
tion  mav  legitimalelv  be  applied  at  microwave  frequencies.  We  will  use  Wright's  expression  |19b(i| 
tor  scattering  from  a  surface  perturbation  when  ka  <  <  I  ( a  is  the  maximum  height  of  the  perturba¬ 
tion): 
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where  the  angle  factors  |.g(  n\~  are  given  by  Valenzuela  |1978),  and  plotted  here  in  Fig.  18  for 
water  at  X-band  frequencies.  For  the  circularly  symmetric  perturbation  in  (8),  the  surface  integral  in 
(9)  becomes  the  Flankel  Transform  of  the  ring  wave  perturbation,  and.  as  shown  in  the  Appendix,  the 
cross  section  given  in  (9)  takes  the  form  given  in  (Alb): 
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Some  typical  values  for  the  ring  wave  parameters  might  he  A=1  and  B  =  .65  as  used  in  Fig.  17. 
V  =  25cm/sec  as  inferred  from  Fig.  16,  a=0.  lent  at  p„  =  lcm.  and  a  =  2k cost/*  =  4  for  X-band 
illumination  at  18  Referring  to  Fig.  18.  we  see  that  for  V-polari/ation.  the  g-factor  is  0.63,  so  put¬ 
ting  all  this  into  ( 19),  we  find  the  ring  wave  cross  section  to  be 

oil)  =0.1 7 sin  tl(K)/  +  7r/4)cm"  (II) 

This  expression  predicts  a  return  that  oscillates  w  ith  a  period  of  about  30  msec  and  with  an  amplitude 
about  6dB  below  the  height  of  the  V-pol  crown  return,  which  almost  exactly  describes  the  small, 
equally  spaced  returns  to  the  right  of  the  main  stalk  return  in  Fig.  6a.  For  horizontal  polarization,  the 

g-factor  given  in  Fig.  18  is  seen  to  be  20dB  lower,  so  the  ring  wave  scattering  for  this  polarization 

would  be  well  below  the  threshold  of  observation. 

VI.  PREDICTION  OF  THE  RAIN-RATE  DEPENDENCE  OF  SPLASH  SCATTER 

The  scattering  models  developed  above  pretty  well  describe  the  scattering  behavior  displayed  in 
the  laboratory  experiments,  but  in  order  to  apply  these  models  to  real  rain  splashes  we  would  have  to 
know  much  more  about  the  characteristics  of  such  splashes  and  the  distributions  of  their  parameters 
than  we  do  at  the  present  time. 

A.  The  Stalk-Height  Distribution: 

We  know'  that  in  natural  rain  the  drops  vary  in  size  according  to  a  drop-size  distribution,  which  is  a 
function  of  rain  rate.  We  have  observed  that  the  splash  dimensions  seem  to  scale,  in  some  way.  with 
the  size  of  the  impacting  drop,  so  we  might  expect  that  the  individual  splash  cross  sections  would  also 

be  distributed  in  the  manner  ot  a  drop-size  distribution.  I  he  stalk  is  clearly  the  most  important 

scattering  feature  in  the  splash,  and  careful  scaling  ot  maximum  stalk  heights  from  video  records  of 
natural  rain  splashes  tor  various  rain  rates  (Cavalcri.  |986|  have  shown  two  cases  in  which  the  distri¬ 
bution  ot  maximum  stalk  heights  looks  verv  much  like  a  l.aws-Parsons  drop-size  distribution.  I  his 
may  be  seen  bv  comparing  the  Cavalcri  stalk-height  distributions  with  the  Laws- Parsons  drop-size  dis¬ 
tributions  | see,  e.g  ,  IMabv.  et  al.  1981 1  for  the  same  ram  rate  R.  as  is  done  in  Fig.  19  for 
R=4mm'hr  and  R  6.4mm/hr.  Such  comparisons  also  yield  an  empirical  estimate  ot  a  scale  (actor 
that  can  he  used  to  convert  drop  sizes  to  splash  dimensions  Note  that  the  peaks  ot  the  two  stalk 
height  distributions  m  Fie.  19  occur  at  h  lk  m  and  h  I  8cm.  respectively ,  while  the  peaks  ot  the 
corresponding  drop-size  distributions  appear  at  I  >  It  and  I  7cm.  ihc  indicated  scale  (actors  ht,/l) 
convert  the  lower  curves  into  the  distributions  sketched  onto  the  upper  curves,  thus  establishing  a 
relationship  between  maximum  sialk  heights  and  the  s|/e  ot  the  drops  producing  them.  While  indi¬ 
vidual  "rains"'  can  have  unique  drop  size  disti ibution  depending  on  meteorological  conditions,  the 
sinnlantv  ot  the  stalk  heniht  and  drop  size  disii ibutions  m  Fig  ll(  suggests  that  the  l.aws-Parsons  dis 
trihution  /i^.il))  might  sometimes  pi  ovule  a  icusonablc  estimate  ot  the  stalk  height  distribution  />(//,,). 
For  such  cases  we  would  vvnte 
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B.  The  Effects  of  Surface  Slopes: 


In  modeling  the  scattering  from  natural  raindrop  splashes  on  a  real  sea  surface,  we  must  make  some 
assumption  about  the  state  of  the  surface.  The  laboratory  measurements  described  earlier  were  made 
with  the  mirror-flat  water  surface  of  an  indoor  tank,  so  the  local  grazing  angle  could  be  put  equal  to 
the  antenna  depression  angle.  But  such  surfaces  arc  quite  unusual  on  open  bodies  of  water,  since 
there  is  usually  a  little  wind  somewhere  on  the  surface  leading  to  some  local  “swell"  or  “chop," 
even  for  “calm"  conditions,  and  in  estuarine  waters  like  the  Chesapeake  Bay,  there  will  be  agitations 
due  to  the  gradients  of  the  tidal  currents  and  to  boating,  both  pleasure  and  commercial.  And  of 
course,  the  splashes  themselves  produce  some  measure  of  “microchop."  For  these  reasons,  the  sur¬ 
face  on  which  the  raindrop  splashes  are  produced  will  generally  have  a  non-vanishing  slope  which 
must  be  taken  into  account  in  defining  the  local  grazing  angle  in  the  scattering  formulas. 

The  grazing  angle  enters  the  expression  for  the  stalk  cross  section  through  the  factors  C  and  F  in  (5). 
In  C  it  appears  in  a  cosine,  whose  small  variation  from  unity  over  the  range  of  angles  of  interest  here 
(2-15°)  will  be  ignored.  The  function  F,  on  the  other  hand,  is  a  sensitive  function  of  grazing  angle, 
through  both  the  kh  sin  yd'  dependence  of  its  trigonometric  functions  and  the  surface  reflection  coeffi¬ 
cient  for  the  vertically  polarized  case.  The  effect  of  surface  slope  on  stalk  cross  section  may  there¬ 
fore  be  calculated  simply  by  replacing  the  |  F  | 2  factor  in  (5)  by  its  average  over  the  surface  slope 
distribution  p(s): 

|  Fv  "  p  =  t°°  |  Fv  "  (*„  +  s)\1p(s)ds  (13) 

where  the  lower  limit  corresponds  to  zero  local  grazing  angle,  below  which  the  stalk  would  be  in 
“shadow."  Although  we  really  do  not  know  what  slope  distribution  to  expect  under  the  conditions 
we  are  discussing,  we  probably  will  not  go  too  far  wrong  by  assuming  it  to  be  Gaussian,  with  rms 
value  ,v„.  Figure  20  illustrates  the  result,  calculated  from  (13),  for  a  grazing  angle  of  2°  and  rms 
slopes  of  0.01,  0.05,  0.10  and  0.15.  Comparison  with  the  corresponding  curves  for  a  2°  grazing 
angle  in  Fig.  9  shows  how  important  the  effects  of  slope  can  be.  particularly  for  horizontal  polariza¬ 
tion  (this  was  also  found  to  be  true  in  the  case  of  sea  scatter  |see  Wetzel,  1987|). 

C.  Calculating  the  Average  Cross  Section: 

The  instantaneous  stalk  cross  section  given  in  (5).  with  the  F- factor  replaced  hv  the  average  (13),  is  a 
function  of  stalk  radius  a.  stalk  height  h  ,  mean  grazing  angle  i/-0  and  radar  frequency  /'.  Both  the 
radius  a  and  the  maximum  stalk  height  ht,  are  observed  to  scale  with  the  drop  diameter  I).  so  it  will 
be  assumed  that  a  is  simply  proportional  to  This  means  that  the  cross  section  in  (5)  may  he  con¬ 
sidered  to  be  a  function  of  maximum  stalk  height  h,,  and  instantaneous  stalk  height  h  . 

os  (h.h„)  ---  \  (’(/».,)■  Ir-  I/-'1  "Un\:  (14) 


where  the  dependence  on  and  I  is  assumed.  There  are  only  two  splash-related  parameters,  h  and 
.  so  if  we  let  pOi R  kllnlh,,  be  the  probability  for  rain  rate  R  that  a  stalk  will  have  a  maximum 
height  between  h„  and  /(,,  r  ,lh,,.  and  be  observed  at  a  height  between  h  and  h  +  dh ,  then  the  aver¬ 
age  cross  section  per  stalk  tor  rain  rate  R  is  written  as 

a  (R)  [  f  a  ill .  It  ,  )/></( .  Ii . :  K  )<lli<lh  , .  1 15) 

•'ll  •>(> 

The  problem,  then,  is  to  find  the  joint  probability  density  />  ( // .  h„  \  K  ). 


.  -v,  w-  w  -v 
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In  discussing  the  laboratory  measurements  of  stalk  scattering  in  section  IV,  we  assumed  that  the  stalk 
rose  and  fell  linearly,  although  there  was  evidence  that  the  stalk  labored  up  to  its  maximum  height, 
and  then  fell  more  quickly.  However,  if  we  again  assume  that  the  rise  and  fall  is  linear,  then  the  dis¬ 
tribution  of  stalk  heights  seen  by  the  radar  will  be  uniform,  so  the  provisional  probability  that  a  stalk 
of  maximum  height  /i„  will  be  observed  at  a  height  h  is  just 


pOi  |/i„)  = 


1  /h„.  h<h„ 
0.  h  >/;„ 


and  the  joint  probability  becomes: 


Pih  ,h,,:R)  —  pill  j  h„  i/h/i,.  ;  R  )  =  ( 1  /h,,  )p(h„:R) .  (17) 

Using  (14)  and  817)  in  (15).  the  average  cross  section  per  stalk  for  rain  rate  R  may  be  written  in  the 
form 


5,1  "</?)  =  J*  Jh,, pih,, :  R )  ./v  H(h„). 


where  all  of  the  scattering  behavior  is  contained  in  the  function 


'v  H(h„)  =  ~  Cv  H(h„)  f '  dh  h:  |  Fv  H(h  )  | 2 

2h,,  11 


By  separating  the  integrand  in  this  way.  it  becomes  possible  to  examine  the  interaction  between  the 
scattering  characteristics  of  the  stalks  and  their  si/e  distributions.  At  the  top  of  Fig.  21  we  show  a 
stalk  height  distribution  pih,,)  of  the  type  found  in  Fig.  19,  while  below  it  the  scattering  functions 
./V  H(h„ )  are  given  for  a  mean  surface  slope  .v„  =  0.1  and  a  stalk  diameter  2a  =  0.1/i„ .  No 
account  has  been  taken  in  the  V-pol  curve  of  the  possibility  of  exciting  the  dielectric-waveguide 
modes  discussed  in  section  IV  B;  the  dashed  line  at  -18  dB  indicates  the  level  associated  with  the 
onset  of  such  a  mode  for  a  stalk  height  corresponding  to  a  drop  size  of  roughly  3mm. 

Taking  the  average  over  the  stalk  height  distribution  according  to  (18).  the  average  stalk  cross  sec¬ 
tions  for  the  6  4mm  hr  rain  rate  become 

(j.1  48  dBnr.  a'1  -  66  dBnr .  R  =  6.4mm/hr,  (20a) 

and  doing  the  same  thing  tor  the  other  rain  rate  in  Fig.  19.  we  obtain 

n'  -  51  dBnr.  a'1  -  75  dBm",  R  -  4.0mm/hr.  (20b) 

This  example  illustrates  what  is  already  ohxious  m  Fig.  21  the  vertically  polarized  return  is  deter¬ 

mined  chielly  by  the  dominant  bulge  in  the  drop  (stalk)  si/e  distribution,  so  is  only  weakly  dependent 
on  the  ram  rate .  while  the  horizontally  polaii/cd  return  is  most  sensitive  to  the  population  of  larger 
drops  (stalks)  in  the  tail  of  the  distribution,  and  is  thus  sensitise  to  both  the  ram  rate  and  the  particu¬ 
lar  shape  ol  the  xi/e  distribution  curve. 

Finally,  in  order  to  convert  mean  cross  sections  per  stalk  into  Normalized  Radar  Cross  Sections 
(NR('S's)  lor  stalk  scattering  in  natural  ram.  we  will  need  the  //;/. i  of  stalk  production  on  the  surface  - 
i.e  .  the  number  ol  stalks  per  square  meter  per  second  produced  on  the  surface,  as  a  function  of  rain 
rate  R  Although  not  every  drop  will  produce  a  stalk,  the  flux  of  raindrops  of  all  sizes  on  the  surface 
is  given  in  Ugai.  et  at.  |  1977|  by  the  empirical  expression 


.  z  '  2/2', V /  •S'/k.-S. 


£ 


/  / 


(21) 


\'n(R)  =  1920  ft"  1  drops /nr  sec 

On  the  other  hand,  Cavaleri  [Id86|  has  measured  the  (lux  of  stalk  production  by  direct  observation, 
and  found  that  while  the  temporal  behavior  of  the  stalk  flux  will  be  different  for  different  “rains.”  a 
reasonable  value  for  rain  rates  of  about  4mm  hr  is  ,V,  —  400  stalks,  nr-sec,  which  is  less  than  a  tenth 
of  the  total  drop  flux  given  by  (21)  for  the  same  R.  Whatever  the  correct  expression  turns  out  to  be, 
the  stalk  flux  must  be  multiplied  by  the  stalk  lifetime  I\  to  get  the  total  number  of  stalks  per  square 
meter  that  will  be  scattering  back  to  the  radar  at  any  instant.  Thus  the  NRCS  for  stalk  scattering 
becomes: 

NRCS:  4=  b,{R)  \\  (/?>/,  nr /nr  (22) 

Unfortunately,  there  are  as  yet  no  measurements  by  which  the  predictions  of  this  theory  could  be 
tested  with  any  confidence.  The  measurements  by  Hansen  several  years  ago  1 1 984) ,  shown  here  in 
Fig.  4,  consisted  of  median  tnot  mean)  backscatter  data  from  the  first  few  minutes  of  a  beginning 
rainfall  Moreover,  it  is  unlikely  that  an  equilibrium  drop-si/e  distribution  would  exist  during  such 
periods,  so  the  kinds  of  statistical  drop  (or  stalk)  size  distributions  we  have  been  discussing  could  not 
even  be  defined,  much  less  parameterized,  by  such  a  number  as  rain  rate  R. 

VII.  SUMMARY  AM)  CONCLUSIONS 

As  a  natural  scatterer.  the  raindrop  splash  has  a  strong  appeal  for  the  scattering  phenomenolo- 
gist.  It  has  three  simple,  cylindrically  symmetric  parts— the  crown,  the  stalk,  and  the  ring  wave. 
These  pans  appear  separately  in  a  well  defined  time  sequence,  and  can  be  treated  as  isolated,  non¬ 
interacting  scattercrs.  following  simple  laws  of  growth  and  decay.  The  crown  and  stalk  resemble  cir¬ 
cular  cylinders  closely  enough  that  well  established  formalisms  may  be  used  to  estimate  their  scatter¬ 
ing  behavior,  and  the  ring  wave  constitutes  one  of  the  few  natural  surface  events  to  which  a  perturba¬ 
tion  theory  may  be  applied  with  any  confidence.  Yet,  laboratory  investigations  of  splash  scattering 
have  disclosed  unexpected  features  in  the  scattering  behavior  that  can  be  explained  only  by  invoking 
more  sophisticated  models  of  the  scattering  process. 

When  applied  to  the  well-defined  and  reproducible  structures  of  a  laboratory  splash,  the  scatter¬ 
ing  models  developed  in  this  report  have  tended  to  describe  the  observed  scattering  behavior  quite 
well.  Although  calibration  ol  the  experimental  returns  was  uncertain,  the  predicted  cross  sections 
appeared  to  be  ol  the  right  order  ol  magnitude.  The  speculative  “dielectric  waveguide”  model  of 
stalk  scattering  lor  vertical  polarization  provided  a  qualitative  explanation  for  the  most  puzzling  aspect 
of  the  laboratory  measurements,  but  while  there  is  some  encouraging  support  for  it  in  recent  measure¬ 
ments.  it  can  be  considered  little  more  than  a  provocative  hypothesis  without  further  experimental 
continuation 

The  major  problems  in  applying  these  models  to  scattering  by  natural  rain  splashes  on  a  real  sea 
surface  lie  m  obtaining  realistic  estimates  of  stalk  heigtht  distributions  and  determining  the  effects  ol 
real  surface  slopes  on  cross  section  predictions  I  -.summing  the  behavior  ol  the  stalk  scattering  func¬ 
tion  against  sample  stalk  height  distributions  leads  to  the  conclusion  that  vertically  polarized  returns 
will  he  onlv  weaklv  dependent  on  ram  rate  while  honzontailv  polarized  returns  will  be  sensitive  to 
both  ram  rate  and  ihe  details  ol  the  height  distribution  curve  Unfortunately,  the  measurements  made 
thus  lar  m  naluiul  tain  l.n  k  the  simultaneous  “mound  truth  essentiallv  drop  (or  stalk  height)  distri 
billions  and  surface  slope  speiita  tuaessarv  lot  meammdul  mtereomparison  between  theory  and 
experiment 


I  I 


The  accumulation  of  a  body  of  experimental  data  that  would  directly  settle  the  outstanding  prob¬ 
lems  in  splash  scattering  on  a  a  natural  water  surface  is  likely  to  be  expensive,  time-consuming,  and 
frustrating.  However,  the  theory  is  really  quite  straight-  forward  and  believable,  so  some  fairly  rea¬ 
sonable  inferences  could  be  made  from  improved  laboratory  studies  involving  smaller  grazing  angles, 
a  wide  assortment  of  drop  sizes,  impact  velocities  closer  to  normal  terminal  velocities,  and  a  surface 
that  could  be  agitated  to  produce  the  varying  slopes  of  a  natural  surface.  Stalk  scattering  via  a 
“dielectric  waveguide"  mode  could  be  investigated  as  a  separate  problem.  Such  information,  coupled 
with  an  expanded  library  of  measured  stalk  height  distributions  and  fluxes  versus  rain  rate,  should 
provide  a  basis  for  an  adequate,  if  yet  somewhat  incomplete,  understanding  of  the  problem  of  rain¬ 
drop  splash  scattering  from  “calm”  sea  surfaces. 

As  mentioned  in  the  Introduction,  this  Report  takes  only  a  first  step  toward  understanding  the 
larger  problem  of  how  rain  affects  scattering  from  the  sea  surface.  Much  work  remains  to  be  done  in 
this  challenging  field. 
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APPENDIX 


1.  Derivation  of  the  Multipath  Factor  F'  11 

In  Fig.  7,  the  scattering  cylinder  is  pictured  as  a  stack  of  discs  of  thickness  dz,  each  scattering 
as  if  it  were  a  slice  of  an  infinite  cylinder  of  radius  a.  Formally,  the  scattered  (magnetic)  field  may 
be  written  in  the  following  integral  form  |Kerr,195l|: 

-ikK  _  . 

Hs(r )  =  -  ik  — -  f  f  Vl,xnxfHr))e  ■  '  dS '  (A  I) 

where  n„  is  the  direction  to  the  observer  and  iixH(r')  is  the  tangential  component  of  the  total  mag¬ 
netic  field  at  the  surface  of  the  cylinder,  and  is  proportional  to  the  incident  field  at  the  point  r '  and  to 
one  or  the  other  of  the  eigenfunction  expansions  in  (1).  (2).  These  two  factors  can  be  separated  by 
writing 


(«,,  v  n  x  H  )  =  In,,  x  n  x  //„  )Hin<  (c  )  (A2) 

where  Hmi.  is  the  combination  of  direct  and  surface-reflected  fields  at  the  height  z,  and  Hlt  contains 
the  eigenfunction  expansion.  For  monostatic  radar,  the  scattered  field  returns  to  the  source  over  the 
same  two  paths  that  delivered  the  incident  field  to  the  disc,  the  two  directions  defining  different  h„'s 
in  (Al).  Although  we  could  go  through  the  formal  manipulations  via  (Al).  it  easier  and  clearer  to 
recognize  that  there  arc  four  possible  ways  a  given  disc  can  interact  with  the  radar:  direct  ray 
incident,  reflected  ray  incident,  scattering  into  direct  and  surface  reflected  rays  in  both  cases.  Let 
Pj  j  denote  the  phases  relative  to  the  reference  plane  RP  in  figure  (Al),  where  1  =  direct  ray. 
2- reflected  ray,  and  RP  is  defined  by  the  grazing  angle  \p  and  maximum  cylinder  height  h. 

The  total  scattered  field  is  obtained  from  two  integrals  of  the  form  of  (Al),  where  only  the 
phase  factors  remain  inside  the  integral  and  all  of  the  remaining  factors  are  contained  in  the  expres¬ 
sion  H(k  .  a  ,  i/\  R )  (verification  is  left  as  an  exercise  for  the  reader): 

H)(IIU  -  Hik.a.t-.R)  jjc'''  ^  Rc1'  +  Re'1’  f  RVr':\dz  (A3) 

R(\p,  f)  is  the  (complex)  reflection  coefficient  for  water  at  frequency  t  and  grazing  angle  i y  and  can 
be  found  in  Saxton  and  Lane  |I955|.  The  phases  are  obtained  from  the  Figure  above,  and  take  the 
form: 
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which  is  used  in  Eqns.  (3)  and  (4)  for  the  metallic  cylinder  cross  sections. 

2.  Reduction  of  Scattering  Integral  for  the  Ring  Wave 

The  scattering  integral  in  Eqn.  (9)  is  given  in  rectangular  coordinates  for  a  plane  wave  incident 
in  the  positive-x  direction: 

l(k  cos  t )  =  j  |  f (.t  ,  \ ,  r  )e ' 2k a,s^ '  dxdy  .  (A6) 

The  ring  wave  is  circularly  symmetric,  so  the  integral  can  be  transformed  into  circular  coordinates 
with  the  substitutions: 

f(.v  ,  v ;  f )  =  f(p:t)  :  x  =  p  cos 6,  y  =  psin0.  dxdy  =  pdpdd  (A7) 

and  written 

l(k  cos \p\  1)  =  \°°  f(p,/)  |f  '  pdp  (AK) 

J()  J() 

00 

=  2tt  f  s"(p;  /)  J„(2k  pcos\P)pdp 
Jo 

which  is  just  2rr  times  the  Hankel  Transform  of  the  ring  wave  perturbation.  Substituting  (8)  into  (A8) 
gives  two  terms  of  the  form: 

(p  -  cm' 

I^(k,\J/\t)  -2ira  p,'/2  j()  e  A  J„(2k  p cos \p) pu2 d p  (A9) 

where  A  corresponds  to  either  A  or  B.  The  ring  wave  is  not  fully  developed  until  it  is  at  least  a  stalk 
height,  or  so,  from  the  splash  center,  so  at  X-band  frequencies  the  argument  of  J„  in  (A9)  will  gen¬ 
erally  be  much  greater  than  unity  over  the  range  of  radii  that  are  important  in  ring  wave  scattering. 
This  means  that  we  can  use  the  asymptotic  approximation  for  the  Bessel  function: 


(A  10) 


J0(ap)  ~  p  ~  sin  (or p  +  rr/4) 

rro: 
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1  |^Mnrp  ♦  Jr/4)  _  ^  /(<» rp  +  ff/4|| 
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=  iv'2tcip,!/-«  1/2  k"r/4  jo  c  A  ci‘"l'dp~c  "/4|(i  r  A  e  ""’dp] 


(All) 


•.*.  •  ■ '  .  *  ,  •  » "  «■»  ."j’<  .‘O,'*  O  ,N  A  Aj'  A  AA,1-.'  .  A  .  .s'.'-  .s'  ■. 


SERIES  I  a — {continued) 


TIME  A 


Separation  of  Droplets 


Dependence  of  Rain  Backscatter  from  Calm  Water 


9  GHZ.  W  SPLIT  SCREEN  VIDEO  DISPLAY 

15  degrees  incidence  STORAGE  POWER  TIM 

SCOPE  (50  ms  div) 

4  mm  drop  - - - — _ _ 


Comparison  of  V-Pol  and  H-Pol  Returns  versus  Drop  Size 


U 


a.  Changing  Field  Configurations  with  Changing  Stalk  Height 


THEORY 


EXPERIMENT 


t  (msec) 


(  a  =  .15mm,  f  =  10  Ghz,  hmax  =  25mm  ) 


b.  Comparison  of  Theory  with  Experiment 


Fig.  14  —  Cross  Section  Variation  with  Stalk  Height  for  Dielectric  Rod  Model 


Fig.  14' 


Theory  and  experiment  for  water  column  of  \arving  height 


Measured 
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